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To study the microscopic electronic and magnetic interactions in the substoichiometric iron chalco- 
genide FeSei-z which is observed to superconduct at x « ~ up to T c =27 K, we use first principles 
methods to study the Se vacancy in this nearly magnetic FeSe system. The vacancy forms a ferri- 
magnetic cluster of eight Fe atoms, which for the ordered x=\ alloy leads to half metallic conduction. 
Similar magnetic clusters are obtained for FeTei-a, and for BaFe2As2 with an As vacancy, although 
neither of these are half metallic. Based on fixed spin density results, we suggest the low energy 
excitations in FeSei-x are antiparamagnon-like with short correlation length. 

PACS numbers: 74.70.-b, 74.25.Jb, 75. 25. +z, 71. 20. Be 



I. INTRODUCTION 

The discovery of superconductivity, [1] now up to 
55 K in the LaFeAsO class of compounds, followed 
by T c of nearly 30 K in the BaFe2As2 class with 
the same active conducting Fe2As2 layer has caused 
tremendous excitement. A competition between su- 
perconductivity and magnetism is very evident, and 
the assortment of phenomena, concepts, and pro- 
posed models is leading to a wide variety of sugges- 
tions that is reminiscent of the heyday of high tem- 
perature superconducting cuprates. There is a seri- 
ous need to identify and address relatively straight- 
forward questions, in addition to broader investiga- 
tions to compare and contrast all the 7?FeAsO ma- 
terials (7?.=rare earth) to identify trends that might 
provide a clue. The system, FeSe, has the same band 
filling as the superconducting oxypnictides, but hav- 
ing only two atoms, is structurally simpler and pro- 
vides more direct questions. Although many samples 
are two phase and are not always fully characterized, 
this system is reported to be non-magnetic and non- 
superconducting at stoichiometry.and is magnetic 
and clearly superconducting 0, 0, m IE at Se sub- 
stoichiometry FeSer-^. 

FeSe can be substoichiometric on either sublat- 
tice and contains two major phases (a and /3). The 
PbO-type a-FeSei-z compound is the one of cur- 
rent interest and has been studied extensively for its 
spintronics-related magnetic properties by Shen and 
coworkers. fA H. lijl [Toj who concluded from observed 
hysteresis a ferromagnetic (FM) state in the nonsto- 
ichiomctric phase, but a nonmagnetic state in the 
stoichiometric phase. Very recently, Hsu et al. re- 
ported superconductivity with T C = 8K, at x = 0.12 
and 0.18. [2| Subsequently, T c has been raised rapidly 
up to 27 K at pressure P = 1.48 GPa,@, H clearly 



putting FeSei_a; in the high T c category with iron 
pnictides having similar band filling. Margadonna 
et al. confirmed T c ~ 14 K at x = 0.08 at am- 
bient pressure. [3| Fang et al. investigated isovalent 
Fe(Sei_j,Tej,)o.82,[E] finding y-dependent T c in the 
range of 8-14 K with maximum at y sa 0.6, dip- 
ping to zero at y = 1. Notably, temperature depen- 
dent susceptibility measurements show an anomaly 
around 100 K, indicating a peculiar and not yet un- 
derstood magnetic instability. 0, Q Hence, the com- 
petition between superconductivity and magnetism 
observed in iron pnictides is clearly extended to these 
iron chalcogenides, and the current picture seems 
to be that superconductivity arises in a phase with 
strong magnetic character. Here we focus on a cru- 
cial feature: Se vacancies are necessary for produc- 
ing the high temperature superconducting state, so 
what is the character of this defect? 



II. RESULTS 

III. STRUCTURE AND CALCULATION 

In the tetragonal phase a-FeSe with PbO (B10) 
structure (space group: PA/nmm, No. 129), Fe and 
Se atoms lie at 2a sites (0,0,0) and at 2c sites (0,^,z), 
respectively. The Fe2Se2 layers have the same struc- 
ture as in LaFeAsO and BaFe2As2- The experimen- 
tal lattice parameters a=3.7693 A and c=5.4847 
A , which are reported recently by Hsu et al. at 
3;=0.12,[2j are used in our calculations. The inter- 
nal parameter z = 0.2372 is optimized by energy 
minimization within the local density approximation 
(LDA). 

Since the superconductivity has been observed 
around x — 0.12 at ambient pressure, we have used 



TABLE I: Optimized structure for a 2x2 supercell with a Se vacancy (space group: P4mm, No. 99), i.e. FegSe7. 
The order of Fe-Fe distances for the relaxed structure is Fel-Fel, Fe2-Fe2, and Fel-Fe2. Fel means Fe atoms near 
a Se vacancy. Note that this structure is optimized in ferromagnetic state. 
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FIG. 1: (Color online) Fermi surfaces (FSs) for nonmagnetic x =0 phase. The FSs, consistent with those shown 
already by Subedi et a/.,[lg| consist of two F-centered hole cylinders, which contain 0.07 and 0.04 holes per a Fc 
respectively, and two compensating Af-centered electron cylinders with more dispersion along the k z direction. 



a 2x2 supercell containing 8 formula units. In this 
supercell, a Se vacancy represents the x = 0.125 
phase, well within the superconducting regime. In 
this phase, this supercell contains two types of Fe 
atoms, one being adjacent to the Se vacancy (Fel) 
and the other farther away (Fe2). As shown in our 
optimized structure given in Table I, the main effect 
of the Se vacancy on the structure is to shift Fel 
atoms toward the vacant site, whereas Fe2 atoms are 



affected little. The changes in interatomic distance 
are -0.14 A for Fel-Fel, +0.07 A for Fel-Fe2, and 
+0.05 A for Fel-Se. 

Two all-electron full-potential codes, FPLO-7 fIil 
Il2| and WIEN2k[l3j based on the augmented plane 
wave+local orbitals (APW+lo) method. [141] have 
been used in these calculations, with consistent re- 
sults. The Brillouin zone was sampled with regu- 
lar dense mesh containing up to 720 irreducible k 
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FIG. 2: (Color online) Total densities of states (DOSs) 
per Fe (both spins), for nonmagnetic x~0 and nonmag- 
netic x = | in the optimized structures, in the regime 
of Fe 3d states. Note that the Fermi level Ef of x = | 
lies on steep side of a sharp peak, promoting a stable 
magnetic state. N(Ef) of x = | at -Ef is 3.07 states/eV 
per one Fe, which is about 3 times larger than that of 
x = 0. 




I I i I i U i I i I i I I 

-3-2-1012 
E-E F (eV) 

FIG. 3: (Color online) Fe atom-projected densities of 
state (PDOSs) of magnetic x — g in the regime of Fe 3d 
states. The vertical line indicates Ef, set to zero, which 
passes through a peak in the minority and nearly bisects 
the gap of ~0.1 eV in the majority. Inset: Blowup PDOS 
in the range of —0.4 to 0.4 eV, clearly showing a gap in 
the majority channel. 



points. Using WIEN2k with the Perdew-Wang LDA 
exchange-correlation functional. [To| the atomic posi- 
tions in the 2x2 supercell with one Se vacancy were 
optimized until forces were smaller than 2 mRy/a.u.. 
For WIEN2k, local orbitals were added to gain flex- 
ibility in dealing with semicore states, Fe 3p and Se 
3e?4s. The basis size was determined by R m tK maa; = 
6. Atomic radii used were 2.21 a.u. for Fe and 1.96 
a.u. for Se. 



IV. RESULTS 

A. Stoichiometric phase. 

All attempts to obtain either FM or antiferromag- 
netic (AFM) states led only to nonmagnetic (NM) 
solution once the Se position is optimized. Even in 
the fixed spin moment calculations, no (meta)stable 
FM state is obtained. The nonmagnetic ground 
state at x = is consistent with experimental 
observations, flfjj Subedi et aZ.[l6[ reported an AFM 
ground state with a very small stabilization energy; 
this result is not in serious conflict since energies and 
moments in these materials are known to be very 
sensitive to structural and computational details. [17| 
We have found such an AF ordering is stable us- 
ing the generalized gradient approximation, but not 
with our LDA approach (both with WIEN2K and 
FPLO). All these results taken together indicate 
FeSe is very near a magnetic critical point. 



As seen in other superconductors containing FeAs 
layers, a transition from tetragonal to a low tem- 
perature orthorhombic structure has been observed 
around T — 70 K.0, H[ As expected from the tiny 
changes in crystal structures, the calculated change 
in electronic structures is slight. We will address 
only the tetragonal phase here. 

The main difference in the band structure of FeSe 
with respect to iron pnictide compounds with 
similar structure occurs along the T — Z line (k z 
direction). FeSe has a similar band structure to 
BaFe9 As9 [l8j . with a flat band (a d orbital lying 
in the xy-pl&ne with a large Fe-Fe hopping integral) 
(at —20 meV) just below Ef, whereas this band lies 
above E F in LaOFeAs[l7j], FeTepJ] and LiFeAs[ll|. 
As might be expected, the Fermi surfaces shown in 
Fig. [T]are less two-dimensional in FeSe than in iron 
pnictide compounds, which have another layer of 
atoms between FePn layers. 

Figure [5] shows the total density of states for x = 
compared with that for x — 0.125. The density of 
states N(Ep) for x = is small, 25% less than the 
value for LaFeAsO and providing no tendency for a 
FM instability. However, there is a van Hove sin- 
gularity at —50 meV, which is absent in the Fe- 
As superconductors. A difference compared to Fe- 
As compounds is the hybridization gap at —3 eV 
that suggests strong Fe-Se hybridization, consistent 
with about 5% smaller Fe-Se distance in this com- 
pound than Fe-As distance in either LaFeAsO or 
BaFe2 As2 . The Se p states also lie somewhat lower 
than the As p states. 
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FIG. 4: (Color online) Blowup band structure of mag- 
netic x = g , showing half-metallicity with total spin mo- 
ment of 0.5 \xb per Fe, near Ef which is set to zero. 



B. Effect of a Se vacancy 

This Se vacancy cannot be treated well by the 
virtual crystal approximation, which we have con- 
firmed by calculations, thus necessitating a super- 
cell treatment of the actual vacancy. The virtual 
crystal approximation replaces both Se atoms and 
the vacancy by a peculiar average entity which ar- 
tificially restores periodicity. Such a model cannot 
bear any relation to the defected system we study. 
For x = 0.125, we find that a strong magnetic state 
centered on the vacancy is stabilized (see below). 
This relaxed magnetic defect, with structure given 
in Table I, gains 32 meV/Fe by structural relaxation. 
Disregarding magnetism, both structures (relaxed 
and unrelaxed) are nearly degenerate. This differ- 
ence reflects important magnetostructural coupling, 
as already found for LaFeAsO,[ljJ and dependence 
of optimized structure on magnetic states is observed 
commonly in the FeAs-based superconductors, [il^, 
[20j which have large calculated spin moment. In 
this relaxed structure, the magnetization energy (en- 
ergy difference between NM and FM) is 133 meV/Fe, 
one-third larger than in the unrelaxed structure and 
reflecting a very strongly magnetic cluster. 

Now we will address the unusual properties of the 
magnetic state. The most interesting point is that 
the two types of Fe ions are aligned antiparallel: Fel 
with 2.14 hb and Fe2 with —1.10 /is in the relaxed 
structure, which reflects a local AFM coupling (as 
opposed to a dclocalizcd spin density wave mecha- 
nism) . This distribution of Fe d states of each spin 
can be observed in the atom-projected densities of 
states given in Fig. [3j The net moment is 0.5 fi,B /Fe, 
and in addition this ordered system is half-metallic 




FIG. 5: (Color online) Total energy differences (solid 
lines, scale on the left axis) and ratio of local moments 
on Fe ions (dashed lines, scale on the right axis) ver- 
sus fixed spin moment, using FPLO and WIEN2k. For 
M psO.75-1 [ib, two distinct states can be found. The 
dashed area is a boundary separating regimes of antipar- 
allel and parallel spin moments (at the boundary, Fe2 has 
nearly zero moment), described by arrows. The symbol 
* at M=0 denotes nonmagnetic state, which has zero 
energy E(0) in this plot. The M—0 state at —65 meV 
has compensating moments on Fel and Fe2 of magnitude 
~1.9 fiB- (For details, see text.) 



(1 /ig/Fe pair), as shown clearly in the band struc- 
ture given in Fig. 01 In the unrelaxed structure each 
Fe ion has smaller moment in magnitude by 0.25 lib, 
although the total moment remains unchanged. The 
Se vacancy leads to creation of antialigned spin mo- 
ments rather than any identifiable charge difference 
between Fe ions. 



C. Fixed spin moment studies 

One may ask: how stable is this antialigned local 
spin state? Fixed spin moment calculations [2lj in 
the FegSe7 compound are used to investigate this 
question. Results can be seen in Fig. [5l which 
shows the expected energy minimum around the 
half- metallic solution with 0.5 /^s/Fe (4 lib per mag- 
netic cluster). Note that the E(M) curve is not 
smooth at a point of half mctallicity[22j where a sub- 
stantial range of "applied field" leads to the same 
unchanging moment. For total moment M=0, two 
solutions are found: the antialigned spin state with 
net zero moment and the simple nonmagnetic state. 
The antialigned state has lower energy by 65 (37 
in FPLO) meV/Fe than the nonmagnetic state, but 
higher energy by 71 meV/Fe than the half metallic 
state. However, the nonmagnetic solution is more 



4 



TABLE II: Optimized structure for a 2x2 supercell with an As vacancy (space group: Pi/mmm, No. 123), i.e. 
Ba4FegAs7, in ferromagnetic state. The lattice parameters [24[ a=3.9625 A and c=12.0168 A were used for the 
optimization. In addition to atoms given below, Ba atoms lie at la site (0,0,0), lc site (5,5,0), and 2/ sites (0,|,0). 
The order of Fe-Fe distances for the relaxed structure is Fel-Fel, Fe2-Fe2, and Fel-Fe2. Fel means Fe atoms near 
an As vacancy. 
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stable than any solution with parallel aligned spins. 
These results indicate that parallel spins are strongly 
antagonistic for this magnetic cluster; the local cou- 
pling is AFM. 

For antialigned spin states, the difference in total 
energy between the two codes used here is associ- 
ated with different local moments on Fe ions. These 
differences between two all-electron, full potential 
codes which usually give equivalent results reaffirm 
the strong sensitivity of the system to small effects. 
At M=0 and 0.25 hb, the moment of Fel obtained in 
WIEN2k is about 15% larger than in FPLO, though 
the ratio of Fe local moments is nearly identical. 
These differences probably reflect the sensitivity of 
FeSei_ x -the softness of its magnetism-to small com- 
putational details rather than representing distinct 
magnetic states. Additionally, changing the total 
moment for the antialigned spin states, the Fel mo- 
ment is insensitive, with only a maximum change of 
10%, whereas the Fe2 moment varies rapidly. From 
such behavior we can conclude that low energy ex- 
citations involve essentially little change in the Fel 
magnetic moment. 

D. Analogies in BaFe2As 2 and FeTe 

To check the robustness of this vacancy induced 
magnetic cluster, we carried out analogous cal- 
culations for FeTei-z and BaFe2As2. A similar 
structural relaxation was performed for FeTeo.875 
as in FeSe,[23| and an As vacancy in BaFe2As2, 
i.e. Ba4FesAs7 given in Table HT1 For BaFe2Asi.75 
and FeTeo.875, the relaxation and magnetic cluster 
are similar although the final states are not half 
metallic, with total moment of 0.7 /is/Fe and 0.42 
fis/Fe, respectively. These magnetic states are fa- 
vored energetically over the nonmagnetic state by 



about 160 meV/Fe for BaFe 2 Asi. 75 and 175 meV/Fe 
for FeTeo.875- 

V. DISCUSSION AND SUMMARY 

Now we consider the broader context. Magnetism 
in superconducting samples, and its possible connec- 
tion to superconductivity, is one of the primary is- 
sues in iron-pnictide superconductivity, and our cal- 
culations establish that Se (or Te, or As) vacancies 
promote strong magnetic clusters surrounding the 
vacancy. Superconductivity occurs only in substoi- 
chiometric samples, and we obtain strong magnetic 
behavior only around Se vacancies. Our fixed spin 
moment results indicate the low energy excitations 
will involve fluctuations in the magnitude of the next 
neighbor Fe spin (relative to the vacancy) , while the 
near neighbor spin remains rigid and antialigned. 
The character of this excitation is antiparamagon- 
like but with short correlation length, a scenario that 
also seems relevant for the iron pnictide supercon- 
ductors. 

In the superconducting FeSeo.88 materials there 
is a magnetic transition characterized by a sharp 
upturn in the susceptibility (apparently also with 
a structural aspect) near 105 K, followed by another 
transition at 75 K where the susceptibility abruptly 
returns to its higher temperature value. @, @ These 
anomalies have not been discussed much yet, but 
the strong magnetic character, and the difference 
in field-cooled and zero-field-cooled susceptibility 
at lower temperature may be reflecting complex 
cluster-glass behavior arising from the immobile but 
interacting magnetic defects that we have studied. 
The appearance of superconductivity in a disordered 
magnetic system such as this provides strong justifi- 
cation for further study of the physics of the FeSei-a; 
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system. 
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